New quantum phases in a one-dimensional Josephson array 
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We examine the phase diagram of an ordered one-dimensional Josephson array of small grains. 
The average grain charge in such a system can be tuned by means of gate voltage. At small grain-to- 
grain conductance, this system is strongly correlated because of the charge discreteness constraint 
(Coulomb blockade). At the gate voltages in the vicinity of the charge degeneracy points, we find 
new phases equivalent to a commensurate charge density wave and to a repulsive Luttinger liquid. 
The existence of these phases can be probed through a special dependence of the Josephson current 
on the gate voltage. 

PACS numbers: 73.20.Fz,03.65.Sq, 05.45. +b 



In the physics of phase transitions, the special class of 
transitions occuring at zero temperature have attracted 
significant attention in the recent years (!]] . At zero tem- 
perature, long-range order may be destroyed by quan- 
tum fluctuations. In some systems, the magnitude of 
these fluctuations is controlled by an external parame- 
ter. Therefore such a system may undergo transitions be- 
tween the phases with qualitatively different behavior of 
long-range correlations. In contrast to the conventional 
finite-temperature transitions, the quantum transition at 
T = can occur even in a one-dimensional system. For 
example, as it was predicted by Bradley and Doniach 
0, a one-dimensional array of superconducting grains 
linked by Josephson junctions may behave as a super- 
conductor or insulator depending on the ratio of charg- 
ing and Josephson energies. Recently a special sample 
design allowing one to vary continuously the strength of 
inter-grain Josephson junctions was developed || and an 
attempt to examine the vicinity of the quantum phase 
transition was undertaken 0] . 

Paper Q and the majority of later theoretical papers 
devoted to the one-dimensional Josephson junction ar- 
rays, assumed that the grains of an array are tuned to 
the point of particle-hole symmetry. In this case, the 
transition from a superconductor to Mott insulator is 
of the Berezinskii-Kosterlitz-Thouless (BKT) type [§. 
Within a broader framework of quantum phase transi- 
tions in a system of interacting bosons it became 
clear that breaking the particle-hole symmetry immedi- 
ately changes the universality class of the transition. 

Much of the physics of interacting bosons can be ap- 
plied directly to the Josephson junction arrays. For an 
array, a boson is replaced by a Cooper pair. Breaking 
the particle-hole symmetry means simply applying some 
voltage V g to a gate which is electrostatically coupled to 
all grains of the array. Breaking the symmetry leads to 
a symmetric in V g shift of the superconductor-insulator 
transition. Recently || some experimental evidence was 
found for such a symmetric shift for two-dimensional sys- 



tems. 

In this paper, we found new, compared to the pre- 
dicted in Q, phases of a one-dimensional Josephson ar- 
ray. These appear due to the effective repulsion between 
Cooper pairs at different sites of the array. In the plane of 
controlling parameters, which arc the Josephson energy 
Ej and voltage V g , these phases occupy the vicinities of 
the charge degeneracy points. At Ej = the effective re- 
pulsion leads to a period-doubling transition. It separates 
the known Mott phase with homogeneously distributed 
Cooper pairs, from a new phase with alternate occupa- 
tion of the grains ("Neel phase"). As Ej increases, the 
regions of Mott and Neel phases separate from each other, 
giving way to another, intermediate phase. This phase is 
equivalent to a Luttinger liquid with repulsion. Further 
increase of the Josephson coupling brings the system to a 
transition into the superconducting state (Luttinger liq- 
uid with attraction). We found the corresponding transi- 
tion line on the phase diagram {Ej, V g }, see Fig. 1. This 
line can be identified with an abrupt change in the behav- 
ior of the Josephson current through a chain containing 
one weaker link. 

We also found the critical exponents for the transition 
line between the Luttinger liquid and Mott phase. 

To start with, we consider a one-dimensional Joseph- 
son array in the absence of Cooper pair repulsion. The 
simplest Hamiltonian that allows for the superconductor- 
insulator phase transition, has the form: 

Hj = -Ej^2cos(<j) l+ i - fa). (1) 

i 

In comparison with , Hamiltonian ([!]) accounts for 

the influence of the gate voltage, eN oc V g (in the ab- 
sence of charge quantization, eAf is the average grain 
charge induced by the gate); Ec is the charging energy 
of a grain. 
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The two parts of the Hamiltonian (yj) do not commute 
with each other, leading to quantum fluctuations of phase 
<j>i of the superconducting order parameter in each grain. 
To start with, we consider the domain of weak Josephson 
coupling, Ej <C EC-, where the fluctuations are strong. 
In the limit Ej — ► 0, the Josephson array is a Mott insu- 
lator at almost any J\f. The exceptions are the discrete 
points J\f = 2n + 1, where a grain with charges 2ne and 
2(n + l)e has the same energies. Josephson tunneling lifts 
this degeneracy. To analyze the properties of the array in 
the vicinity of the degeneracy points, we can project the 
Hamiltonian (Q) on the subspace of states with charges of 
each grain confined to the values 2ne and 2(n+ l)e only. 
In the lowest order in Ej the projected Hamiltonian has 
the form: 

i 

= * Ek+^+^ix)* (2) 

i 

where a z , a + , o~ are the Pauli matrices, and the "mag- 
netic field" h = (TV — 2n — l)/2 allows for tuning the 
system to a degeneracy point by means of the gate volt- 
age. The Hamiltonian ^ can be converted to the con- 
ventional tight-binding Hamiltonian for free fcrmions by 
means of the Jordan- Wigner transformation, and then 
diagonalized Q in the plane wave representation: 

Hq = E £ fca[a/c, £k = Ej cos k - E c h. (3) 
fc 

The fermion band is completely full or totally empty, 
if \h\ > Ej/Eq- This condition defines the location of 
the Mott phase in the plane of dimensionless parameters 
{W, Ej/Ec}. The Mott gap, e M = \E c \h\-Ej\, co- 
incides with the lowest possible excitation energy, and 
vanishes at the phase boundary. Clearly, the phase di- 
agram is periodic in Af with period 2. This allows us 
to consider further only a "strip" of the phase diagram, 
< J\f < 2 which corresponds to n = in the above 
formulas. 

The lowest order in Ej we used is sufficient to find 
the characteristics of the Mott phase. The complemen- 
tary phase in this approximation is equivalent to a non- 
interacting Fermi gas. It is known however, that even a 
weak interaction is relevant in one dimension, transform- 
ing the gas into a Luttinger liquid. There are two com- 
peting mechanisms of interaction between the fermions. 

The inter-grain electrostatic interaction leads to an ef- 
fective repulsion between the Cooper pairs. We will con- 
centrate on the case of a relatively weak inter-grain inter- 
action, E z <C Ec, which allows us to use Hamiltonian (|^) 
as a starting point. In a realistic system with a gate, the 
interaction between distant grains is screened out, and 
only the nearest-neighbor term should be retained: 



H z = J2E z *t<T? +1 . (4) 

i 

On the other hand, the higher order expansion in Ej/Eq 
leads to an effective attraction between the Jordan- 
Wigner fermions. To find the corresponding correction 
to the low-energy Hamiltonian (||), we should take into 
account the virtual states with energies exceeding Ec ■ To 
do this, we use the perturbation theory. In the second- 
order perturbation theory, the effective Hamiltonian re- 
duced to the subspace of charges 0, 2, takes the form: 

H 2 = PHj(l - P)—±—HjP. (5) 

e — He 

Here P = \\ i Pi, an d Pi is the projection operator on 
the subspace of states with charges and 2 for each 
grain; energy e belongs to the reduced energy band, and 
may be neglected in comparison with Ec- To implement 
the projection procedure in (|5|), we use the identities, 
P i exp(±0 i )(l - Pi)exp{T<Pi)Pi = (IT of)/2, that can 
be easily checked. Thus, the result of projection can be 
written in terms of Pauli matrices, 

= -if- E + + -ti-r + i) ■ (6) 

The full Hamiltonian + H z + H 2 enables us to in- 
vestigate the phase diagram in the vicinity of the charge 
degeneracy line (h <C 1). 

We can neglect the interaction H 2 , if Ej/EcE z <C 1; 
in this limit, our system with Hamiltonian + H z 
is equivalent to a one-dimensional anisotropic Heisen- 
berg model in a magnetic field. It is known JTo| that 
at Ej < AE Z and relatively small field \h\ < h%, the an- 
tiferromagnetic Izing interaction dominates the physics 
of the anisotropic Hcscnbcrg model, whereas in a large 
field, \h\ > ft,2, the chain is in the ferromagnetic state. 
Returning to the language of the initial chain of grains, 
the ferromagnetic phase corresponds to Mott insulator. 
The low-field Neel phase is a commensurate charge den- 
sity wave state with period 2 (on average, every second 
grain of the chain is occupied by a Cooper pair). In both 
states, there is a gap for excitations. The region in be- 
tween, h\ < \h\ < h 2 , is occupied by Luttinger liquid. 
The width of the Luttinger liquid region depends on the 
value of Ej. In the limiting case Ej = which corre- 
sponds to the Izing model, this width shrinks to zero. 
At small Ej, the dependence of both hi and h 2 is lin- 
ear in Ej with coefficients of opposite sign jl0| , and the 
width of the Luttinger liquid region is proportional to Ej. 
The field hi turns zero for an isotropic antiferromagnet, 
Ej = AEc- At larger Ej, the Neel state is absent, and 
the liquid phase exists even at h = 0, see Fig. 1. 

With the further increase of Ej, the interaction H 2 
becomes important. As we will see, it alters the sign of 
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interaction in the Luttinger liquid. To study the inter- 
action, we transform the spin Hamiltonian H z + Hi to 
the Jordan- Wigner fermion representation. In addition 
to quadratic in a, terms, this Hamiltonian gives rise 
also to two-body interaction Vj n t • The significant part of 
(weak) interaction V ln t that defines the properties of Lut- 
tinger liquid, corresponds to the interaction between the 
left- and right-movers. Near the one-dimensional Fermi 
"surface" k = ±kp this part can be written as: 

Vint = EjX ^ sin2 k F a} ki a} k2 a k2+q a kl ^ q , (7) 
fci ,fc 2 ,<? 

A= f 8Ez 7E A 

V Ej 2E C ) ' 

Here k\ f» kp, k% ~ —kp, and q <C kp; the Fermi 
wavevector varies with the gate voltage. 

Interaction ((?]) vanishes at the phase boundary of the 
Luttinger liquid with Mott phase (kp — > or kp — > tt), 
regardless the relations between parameters E z , Ej, and 
Ec- Therefore, we can use the free- fermion model to 
study the critical behavior of our system near the phase 
boundary. The divergent length scale, £(h) ~ 1/kp, char- 
acterizes the distance between free charges in the con- 
ducting phase. Using Eq. (|3|), we find £ oc S^ 1 ^ 2 , so that 
the critical exponent v equals 1/2 (here 8 = hi — is 
the distance to the phase boundary). The characterisitic 
frequency scale fi is given simply by the Fermi energy. 
The corresponding exponent z, defined by the relation 
Q ~ 5 ZU , is z = 2. Finally the superfluid density ||, 
p s oc 5^ , should be identified with the density of fermions 
p s ~ 5 1 / 2 , which determines £ — 1/2. 

Vanishing of the interaction at the phase boundary 
with Mott phase is a property far more general than the 
result of perturbation theory ([?]) , and it reflects the fact 
that we deal with spinless fermions. In the higher orders 
in Ej/Eq the dependence of hi on Ej becomes non- 
linear, and the spectrum of fermions changes, but the 
critical exponents we found are not affected as long as 
TV 7^ 0, 2. The lines TV = 0, 2 are the lines of particle-hole 
symmetry. Because of this symmetry, the system belongs 
to a different universality class || , and the transition to 
the Mott phase (which occurs at Ej/Eq ~ 4/-7T 2 ) is of 
the BKT type. 

Now we discuss the properties of the Josephson junc- 
tions chain in the Luttinger liquid state, which low- 
energy properties are described by a quadratic La- 
grangian y|§, 

C=^-[v(d x 4>) 2 +v- 1 (d T cf>) 2 ]. (8) 

Here v is the velocity of the acoustic excitations in the 
system (plasmons in our case), and g determines the 
long-range behavior of all correlation functions. The phe- 
nomenological constant g can be calculated in terms of 
the microscopic parameters in some limiting cases. 




FIG. 1. Zero-temperature phases of the one-dimensional 
array of Josephson junctions controlled by the gate voltage 
(V g ) and Josephson energy (Ej). In addition to the known p] 
Mott phase M and superconducting phase S (which is equiv- 
alent to the attractive Luttinger liquid), we find the commen- 
surate charge density wave state N, and repulsive Luttinger 
liquid phase RL. The Luttinger liquid parameter g equals 1 
at the RL-S boundary, and also at the boundary of the Mott 
phase; g — 1/4 at RL-N boundary. The special points on 
the phase transition line where the charge degeneracy occurs 
(a), and where the particle-hole symmetry is observed (6), are 
multicritical points with g = 1/2 and g = 2 respectively. 

Deep in the superconducting region, Ej ^> Ec, one 
can compare Eq. (j^) with the quadratic Lagrangian de- 
rived by expansion of (|l|) in <ft, to find g « tt(Ej /Ec) 1 ^ 2 ■ 
In the opposite limit of a relatively weak Josephson cou- 
pling, E z Ej <C Ec, constant g can be found with the 
help of (@), g w 1 - (A/2tt) sin kp. In addition to these 
limiting expressions for g, we know its value at the phase 
boundaries. We already mentioned that at the boundary 
with the Mott phase the interaction is vanishing, which 
means o = 1. The exceptions are the special multicritical 
points H , where the boundary intersects with a particle- 
hole symmetry line (point b in Fig. 1). If point b is crossed 
along the symmetry line, the transition is of the BKT 
type with g = 2. In the vicinity of these points, g rapidly 
varies from g = 1 to g = 2. With the help of Haldane's 
theory p3[ for anizotropic s = 1/2 Hezenberg model in a 
magnetic field, we find g = 1/4 for the boundary between 
the Luttinger liquid and Neel phases. The exceptions are 
the multicritical points (a in Fig. 1) of charge degeneracy, 
h = 0, where g — 1/2. Similarly to the lines of particle- 
hole symmetry, the phase transition at h = is of the 
BKT type, see, e.g., @. 

Different phases on the diagram Fig. 1 can be distin- 
guished by the dependence of the critical Josephson cur- 
rent on the length of the system L. Clearly, the gap in 
the excitation spectrum existing in the insulating phases 
leads to an exponential suppression of the critical current 
with L] the corresponding correlation length diverges 
(with exponent v) at the transition line. 

To characterize the Luttinger liquid state, we consider 
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a chain of junctions connecting two massive supercon- 
ducting leads, and containing just one especially weak 
link, E w <C Ej. Deep in the superconducting state 
(g 3> 1) we can neglect with the quantum fluctuations. 
The energy of the system consisting of L junctions of the 
nominal strength Ej and one weak link equals: 

m = Ej (<t>-M 2 _ EwCOS(/)w (9) 

Here <f> is the phase difference applied to the leads, and <f> w 
is the phase difference accross the weak link; the energy 
(^) should be minimized with respect to 4> w . Quantum 
fluctuations of phases 4>i at finite g lead to zero-point mo- 
tion of 4> w around some average value <j> w , which depends 
on the external phase cf>. Averaging over the fluctuations 
6cj> w — 4> w — <p w results in the energy functional of the 
form (||) with E w renormalized by a proper Debye- Waller 
factor, and <f> w replaced by <f> w , see, e.g., jl5|. For a short 
array, the renormalized Josephson coupling constant is 
E^ ~ E w (l /L) 1 / 9 . Energy ^ corresponds to two in- 
ductors in series. Therefore, the response of the system 
to the external phase is dominated by the weaker of two 
elements. If g > 1 (weaker quantum fluctuations), the 
renormalized energy £!?? decreases with L slower than 
the first term in Eq. (||). Consequently, at L exceed- 
ing the crossover length, L* ~ (Ej / E w ) 9 ^ 9 ^ x \ the en- 
ergy (||) is dominated by the first term with 4> w = 2nn. 
The Josephson current cx dE/d(f> has a sawtooth de- 
pendence on phase, with the amplitude proportional to 
Ej/L. Thus, the effect of the weak link is "healed" over 
the distance L* . 

If g < 1, the weak link energy E^f falls off faster 
than 1/L, and this link remains the weakest of the two 
inductors connected in series at any L; consequently, 
4> w = 4>. Therefore, "healing" does not occur, and at 
any L the chain behaves as a single Josephson junction 
of strength E%. , . The corresponding current-phase rela- 
tion is I ~ EyjL 1 / 9 sin 4>. 

The line g = 1 on the phase diagram may be viewed 
as a phase transition [[l6f between the superconducting 
and insulating phases in the following sense. In the limit 
L — > oo, the inductance C of an array with a single 
weak link is insensitive to the weak link if g > 1 (su- 
perconducting state), and C oc L/Ej. In the insulating 
state, the scaling of inductance with L is (/-dependent, 
C cx L 1 ^ 9 /E w , with exponent l/g>l @. 

To facilitate the quantitative analysis of the phase di- 
agram, we considered the limit A ^> Eq 3> E z . Eazing 
the first, and lifting the second of these two restrictions 
does not alter the diagram qialitatively. Moreover, in a 
real system with large inter-grain capacitance B the en- 
ergies Eq and E z are of the same order, and the domains 
of charge density wave and of the repulsive Luttinger liq- 
uid are not small. The phase diagram changes drastically 
only if the superconducting gap A becomes smaller than 
some critical value, at which the grains may carry single 



electrons in the ground state. The charge density wave 
state is destoyed then, and a phase appears with an odd 
charge on each grain. For this phase, the Josephson cur- 
rent scales exponentially with L, but there is no gap in 
the excitation spectrum due to the dense quasiparticle 
spectrum in each superconducting grain. 

We assumed a perfect order in the system, which 
means there are no spatial variations of the parameters 
of the Hamiltonian (Q) . The validity of the above analy- 
sis requires that the Josephson junctions chain is shorter 
than the mean free path for a Cooper pair propagating 
along the chain. 

In conclusion, we examined the zero-temperature 
phase diagram of a one-dimensional chain of Josephson 
junctions. We have shown that accounting for the re- 
pulsion between Cooper pairs occupying different grains 
leads to new phases, which are equivalent to the phases 
of one-dimensional spinless fermions with repulsion on a 
lattice. We have studied in detail the boundaries of the 
quantum phases. 
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